High-order harmonic generation and multi-photon ionization of ethylene in laser fields 
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Applying time-dependent local density approximation (TDLDA), we study the high-order har- 
monic generation (HHG) of ethylene subjected to the one-color (w = 2.72 eV) and the two-color 
(coi = 2.72 eV and u>2 = 5.44 eV) ultrashort intense laser pulses. The HHG spectrum of ethylene 
in the one-color laser field shows the obvious plateaus and odd order harmonics are produced while 
the two-color laser field can result in the breaking of the symmetry and generation of the even order 
harmonic. The ionization probabilities are obtained showing the increase of the ionization proba- 
bility of higher charge state by the two-color laser field. The temporal structures of HHG spectrum 
of ethylene is explored by means of the time-frequency analysis showing new insights of the HHG 
mechanisms in the one-color and the two-color laser fields. 

PACS numbers: 42.65, 33.80, 34.80 



I. INTRODUCTION 



The interaction between strong laser and atoms and 
molecules is a hot topic. A series of nonlinear phenom- 
ena appear when atoms and molecules are subjected to 
intense laser pulse p], 0, 0, 13] ■ High-order harmonic gen- 
eration (HHG) is one of the most studied effects among 
these nonlinear phenomena with the results that the co- 
herent emitted light can be harnessed to produce new co- 
herent trains of extremely short pulses on the time scale 
of electron motion [H, [|| . A simple model Q , based on 
the idea of tunnelling ionization of an electron in an atom 
allows for interpreting the harmonic generation process 
in terms of recolliding electron trajectories and leads to 
a universal cutoff law for the maximum in photon energy 
in harmonic generation. The extra degrees of freedom 
and nonspherical symmetry of molecules lead to a fur- 
ther new phenomena such as molecular high-order har- 
monic generation (MHOHG) JEMS, charge resonance 
enhanced ionization (CREI) [fj, and bond softening 
[TTI . EH. The theoretical understanding of these mech- 
anisms needs to solute the time-dependent Schrodingcr 
equation for all electrons and all nuclear degrees of free- 
dom. However, this kind of numerical solutions only exist 
for smallest systems, such as atoms [H|,EI|], H2 and H^ 
[HI, E3, EH, El, EH ■ For large systems the computation 
is quite expensive due to the number of degrees of free- 
dom. In contrast, the well tested time dependent Den- 
sity Functional Theory at the level of the time-dependent 
local-density approximation (TDLDA) [2(j serves as a 
powerful tool to stud y th e electron dynamics of multi- 
electron systems [2ll l22l EH- Recently, Madsen et al 



[24| studied MHOHG of polyatomic molecules using a 
quantum-mechanical three-step model. 

The application of two-color laser field for HHG is a 
fascinating topic of research science. One can control 
the formation of harmonic spectrum and yields of ion 
or molecular fragments by changing the relative phase 
between the fundamental frequency and its second har- 
monic fields in the two-color laser field. At present, in 
experiment and theory, two-color laser field has been ap- 
plied to atomsEIJMjp, molecules [H, Eil H EH, E3 
and clusters [33[~34, 35|. However, much less has been 
done for more complicated molecules. In this paper we 
apply TDLDA, augmented with an average-density self- 
interaction correction (ADSIC) [36] to investigate theo- 
retically the ionization and MHOHG of ethylene in one- 
color and two-color laser fields. 



II. THEORY 

In this section, we represent the real-time method in 
TDLDA. The molecule is described as a system com- 
posed of valence electrons and ions. The interaction be- 
tween ions and electrons is described by means of a norm- 
conserving pscudopotentials. 

Valence electrons are treated by TDLDA, augmented 
with an average-density self-interaction correction (AD- 
SIC) [36]. They are represented by single-particle or- 
bitals (r, t) satisfying the time-dependent Kohn-Sham 



(TDKS) equation [3J 
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V e f / is Kohn-Sham effective potential composed of four 
parts, 

V eff (T,t)) =V ion (r,t) + V ex t(r,t) 

+V H [n](r,t) + V xc [n](T,t), (2) 

where Vi on = X)/^ps( r — ^r) * s i° n i c background po- 
tential, V ex t is external potential, Vh stands for a time- 
dependent Hartree part and the final part is exchange- 
correlation (xc) potential. The electron density is given 
by 

n(r,t)) = £>,-(M)| 2 , (3) 
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and the Hartree potential V# [n](r, t) is defined as 

V H [n](r,t)=Jd*r>^±. (4) 

The xc potential 14 c [n](r,i) is a functional of the time- 
dependent density and has to be approximated in prac- 
tice. The simplest choice consists in the TDLDA, defined 
as 

V x T c DLDA [n](r,t) = de h x r(n)/dn\ n=n{r , th (5) 

where e^° m (n) is the xc energy density of the homoge- 
neous electron gas. For e^° m we use the parametrization 
of Perdew and Zunger [38[ . The form of pseudopotential 
for covalent molecule is taken from [3^ | including nonlo- 
cal part. 

The ground state wavefunctions are determined by the 
damped gradient method |4(j. The TDLDA equations 
are solved numerically by time-splitting technique [4lj . 
For the nonlocal part contained in Hamiltonian, we deal 
with it in an additional propagator and treated it with 
a third-order Taylor expansion of the exponential [iol ]. 
The absorbing boundary condition is employed to avoid 
periodic reflecting electrons [42j ■ 

For the laser field, neglecting the magnetic field com- 
ponent can be written as 

E(t) = E f(t)[cos(Lu{t)) + Acos(2uj(t) + 0)] (6) 

where Eq cx v7, I denoting the laser intensity, u> is the 
laser frequency and fit) is the pulse profile. A is the 
electric field strength ration between the two frequency 
and 4> is the relative phase. In this paper, A is chosen 
and 0.5 respectively for the one-color and the two-color 
laser field and = 0. 

The laser induced electron dipole moment can be ob- 
tained by 

D(t) = J rn(r,t)d 3 r (7) 

The Fourier transform of D(t) gives the MHOHG power 
spectrum \Dp(uj)\ 2 . 



The harmonics intensity as a function of harmonic fre- 
quency lo and emission time (3 can be obtained by the 
time-frequency analysis [43| 

D g {uj,P) = | J D(t)e l ^e- {t - p)2/2a2 dt\ 2 (8) 

where the window function width a is chosen as one-tenth 
of the laser optical period. 

The number of escaped electrons is defined as 

N esc = N t = - [ d 3 rn(r,t) (9) 
Jv 

where V is a volume surrounding molecule. A detailed 
link with experiments is the probabilities P k (t) of finding 
the excited molecules in one of the possible charge states 
k to which they can ionize. The formula can be obtained 
from [42| . 

III. RESULTS AND DISCUSSION 

The ethylene is the simplest organic n system holding 
D2h symmetry. In our calculation, there are 12 valence 
electrons. It is in x-y plane with the center of mass at the 
origin. The laser polarization is along x direction which 
is parallel to the axis of CC double bond. The laser 
time profile is chosen the Ramp envelope. The ramp- 
on and ramp-off time are both 8 fs and the duration is 
30 fs. For the one-color laser field, I = 10 14 W/cm 2 
and oj = 2.72 eV and for the two-color laser field, I\ = 
10 14 W/cm 2 , I 2 = Ji/4 and uj x = 2.72 eV, lj 2 = 5.44 eV. 
It should be noted that the two-color laser field is non- 
symmetric, which is the same as that in [3l[, the negative 
and positive amplitudes of the electric field strength are 
not equal. 

A. Ionization properties 

Fig. [T] shows the time evolutions of the number of es- 
caped electrons of ethylene in the one-color and the two- 
color laser fields. One can find that the electron emission 
takes the same pattern in both cases. The ionizations 
start at around 5 fs with a constant speed in both cases, 
while it is obvious that electrons start to escape a little 
earlier in the two-color laser field. Then in both cases 
electron emissions become more and more slow. Finally, 
at around 42 fs they are saturated at 1.18 in the two-color 
case and 0.42 in the one-color case, which is 4 fs earlier 
than the laser pulses are switched off. Moreover, it is 
obvious that the ionization is enhanced by the two-color 
laser field. 

Fig. [5] represents ionization probabilities of various 
charge states for ethylene in one-color and two-color laser 
fields. From Fig. [Ha) we can see that in the first 5 fs, 
there is mainly neutral ethylene in the one-color laser 
field and probabilities of _p( 1+ ) and p( 2+ ' arc almost zero. 
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FIG. 1: Time evolutions of the number of escaped electrons of 
ethylene in the one-color and the two-color laser fields. Solid 
line: one-color (u = 2.72 eV) laser field with / = 10 14 W/cm 2 . 
Dashed line: two-color laser field with Wi = 2.72 eV, Ji = 
10 14 W/cm 2 and u>2 = 5.44 eV, h = ii/4, the relative phase 
<?!> = 0. 
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FIG. 2: The ionization probabilities of ethylene in the one- 
color laser field (a) and two-color laser field (b) for the same 
cases as in Fig. [T] 



Then from 5 fs to 42 fs, the probability of neutral ethylene 
decreases more and more slowly and reaches a saturated 
value 0.62 before the laser is switched off. During this 
period, p( 1+ ) and p( 2 +) increase more and more slowly 
and they are also saturated at 42 fs. Their final values are 
0.3 and 0.05 respectively. The probability of p( 3+ ) is so 
small that it is not visible in the figure. Furthermore, one 
can find that in the course of the one-color laser pulse, 
neutral ethylene predominates the main probability. 

For the ionization probabilities of ethylene in the two- 
color laser field, as shown in Fig.[2]Jb), it is obvious that 
the time evolution pattern is different from that in the 
one-color case. In the first 5 fs, there is mainly neutral 
ethylene in the the two-color laser field. From 5 fs to 
42 fs, p(°> decreases more and more slowly and reaches 
a saturated value 0.28. This value is smaller than that 
in Fig. E{a). It is noteworthy that from 5 fs to 10 fs, 
p(o) d r0 p S quickly from 1 to 0.62 in the two-color case 
while it takes 37 fs for p(°) to drop from 1 to 0.62 in 
the one-color case. In Fig. EJb), p( 1+ ) starts to increase 
quickly at around 5 fs and it exceeds P 1 - -* at around 23 fs. 
Finally, it is saturated at 0.4. For P^ 2+ \ it starts to 



increase at around 7 fs, which is a litter later than P^ 1+ \ 
but earlier than P ( - 3+ - ) . The increase trends of p( 2 +) and 
p(3+) arc q m te similar. Both of them increase slowly 
and reach saturated values 0.22 and 0.06 respectively at 
around 42 fs. We can also find that after 23 fs, p( 1+ ) 
dominates mainly the charge state. 

From the above discussion we can say that the prob- 
ability of higher charge state is larger for the two-color 
than that for the one-color laser field. This is related to 
the fact that the local maximum of amplitude strength 
for the two-color laser field is larger than for the one- 
color case. This is consistent with the result from the 
interaction between atoms and laser fields [44| . 

B. Analysis of MHOHG 
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FIG. 3: (a) and (b): the time evolution of dipole moments 
of ethylene along x, y and z directions in the one-color and 
the two-color laser fields, (c) and (d): MHOHG spectrum of 
ethylene in the one-color and two-color laser fields. The laser 
parameters are the same as those in Fig. □ 

Figs. G2Ja) and (b) exhibit the time evolution of dipole 
moments along x, y and z directions of ethylene in the 
one-color (w) and the two-color (lu + 2uj) laser field re- 
spectively. One can find that in both cases the dipole 
moments along x direction are distinct while those along 
y and z directions are very small. This is due to the laser 
polarization is along x direction. However, D x displays 
the different behaviors in two cases. In Fig. [3{a), D x ex- 
hibits the relative symmetry while in Fig. [3Jb) it is not 
symmetric because of the non-symmetric property of the 
electric field strength of the two-color laser field. 

The MHOHG spectra calculated by our TDLDA in 
the one-color and the two-color laser fields are shown 
in Figs. Etc) and Fig. Eld). It is well known that the 
chosen laser intensities and frequencies are certain essen- 
tial physical parameters which allow for quasistatic in- 
terpretations of strong field laser-atom processes 0, HH . 
In this paper, ponderomotive energy U p = el/4muj 2 for 
u = 2.72 eV and 2uj = 5.44 eV are 1.93 eV and 0.12 eV 
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respectively. The Keldysh parameter 7 separating mul- 
tiphoton and tunnelling ionization regimes for ethylene 
ionization potential I p = 11.5 cV is 7 = y // I p /2U p = 1.7. 
Thus the present parameters situate our calculations 
above the tunnelling ionization regime. 

According to the atomic HG cut-off law, E max = 
I p + 3.17U P , the cutoff MHOHG of ethylene in the one- 
color laser field is at 6. From the MHOHG spectrum in 
Fig. [3jc) we can find a clearly cut-off at the 9th order. 
Although this value is not exactly the same as the value 
predicted by the atomic HG cut-off law, it is quite near 
and we attribute the difference to the fact that ethylene 
is polyatomic molecule while the cut-off law is for atoms. 
One can also find in Fig. [3{c) that the MHOHG spec- 
trum of ethylene consists of a series of peaks, first de- 
creasing in amplitude and then reaching a plateau. The 
plateaus is rather short from the 5th order harmonic to 
the 9th order followed by continuously decrease. Thus 
the MHOHG spectrum of ethylene behaves basically like 
an atom with rccollision of the electron with the com- 
pact molecular ion being the principal HG mechanism. 
Furthermore, one can find in Fig. [3Jc) that ethylene gen- 
erates odd harmonics. 

Compare the harmonic spectrum in Fig. [3jd) to that 
in Fig. 02c), we can see that the shapes of MHOHG spec- 
trum are different. In Fig. [2Id), the plateaus is not so 
obvious and the decrease lasts longer. This is due to 
the interference effect of the lw and 2uj components of 
the two-color laser field, ft can also be found that even 
harmonics are well produced in the two-color case. One 
explanation of this is that even harmonic results from 
the sum of an odd number of 2uj photon plus two lto 
photons. The other interpretation is that even harmonic 
generation requires broken reflection symmetry because 
only this allows that even harmonic generation trans- 
forms a squared dipole excitation into one dipole signal, 
i.e., D 2 — >D. That transition cannot be mediated by a 
reflection symmetric system because parity is then con- 
served, but D has negative parity while D 2 has positive 
parity. Ethylene is symmetric and a two-color laser field 
breaks the symmetry so that a series of even harmonics 
are produced. 

To investigate the detailed temporal structure of the 
MHOHG spectrum of ethylene, we perform the Gabor 
transform of the dipole moments in the one-color and 
the two-color laser field respectively. The laser parame- 
ters are the same as those in Fig. [T] Here we consider two 
typical harmonics which are both in the plateau, one is 
odd harmonic (5th) and the other one is even harmonic 
(6th), as shown in Fig. 2J It should be noted that the 
optical cycle refers to the one-color laser frequency. In 
Fig. HJa) we can find that in the one-color laser field, the 
time profile of the 5th harmonic shows a relative smooth 
function of the driving laser pulse which indicates that 
the multiphoton mechanism dominates this lower har- 
monic regime. However, for the two-color case, the time 
profile of the 5th harmonic is stronger and it exhibits well 
the one burst within each optical cycle from the 5th opti- 
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FIG. 4: Time profiles of the harmonics at about the 5th one 
and the 6th one in the one-color laser field (solid line) and the 
two-color laser field (dotted line). The laser parameters are 
the same as those in Fig. [TJ 



cal cycle to the 25th optical cycle. This can be attributed 
to the recollision of the electronic wave packet with the 
ionic cores. Due to the first five and the final five opti- 
cal cycles are ramp-on and ramp-off times of laser pulses, 
there are only two and four bursts in these two parts of 
times. For the 6th harmonic, as shown in Fig.QJb), it is 
obvious that in the one-color laser field, the time profile of 
the 6th harmonic is very weak, whereas in the two-color 
laser field, it is almost two orders of magnitude stronger. 
Again, it shows well the one burst within each optical 
cycle from the 5th optical cycle to the 25th optical cycle. 



IV. CONCLUSIONS 



In this paper, we have simulated the MHOHG spectra 
and ionization of ethylene induced by the one-color and 
the two-color laser fields in the multiphoton regimes with 
TDLDA. We find that the ionization and the ionization 
probability of higher charge state are enhanced by the 
two-color laser field. It is shown that the MHOHG spec- 
trum of ethylene in the one-color laser field exhibits the 
typical atom HHG spectrum and odd order harmonics 
are produced. The two-color laser field can result in the 
breaking of the symmetry and generate the even order 
harmonics. Furthermore, the detailed temporal structure 
of MHOHG spectrum of ethylene is obtained by means 
of the time- frequency transform providing new insights of 
the MHOHG mechanisms in the one-color and two-color 
laser fields. 
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